Method for the ground-based monitoring of EWF-type anomalies in a positioning satellite signal

ABSTRACT

The method of the invention, which is applicable, in particular, to making positioning signals from GPS and SBAS satellites reliable, consists in performing an instantaneous statistical analysis (of around 1 mn duration) of the correlation peak from the satellite signal receiver and in comparing it to a long-term (several hours) statistical analysis by choosing 5 points which are characteristic of the peak (prompt, ±0.044 and ±0.088). If the result of the comparison exceeds a given threshold, the corresponding signals are rejected.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to a method for the ground-basedmonitoring of EWF-type anomalies in a positioning satellite signal.

DESCRIPTION OF THE RELATED ART

2. Discussion of the Background

Currently, aircraft landing guidance systems are of the ILS type, but,for economic reasons, various States are looking to replace them withGLS guidance systems that use information supplied by satellitenetworks, in particular GPS. In aircraft, the on-board part of theselanding approach guidance systems is of the MMR type that combines ILSsystems with GLS and MLS systems.

GLS systems would be the most economical, especially in view of the factthat GPS positioning can supply aircraft with the information needed fortheir navigation. The performance levels required for navigation incruise flight allow the GPS system to be used in an autonomous manneroffering an accuracy of around 20 to 30 meters and an integritysufficient for the requirements. In the landing approach phase, wherethe required vertical precision is around 2 to 3 meters, the system isused in combination with a complementary ground-based means thatprovides the information necessary to improve the precision and makesthe mechanisms available that allow the positioning integrity to beguaranteed.

Numerous studies have been undertaken to determine the sources of thepositioning errors affecting GPS measurements. The intentionaldegradation of the measurements carried out at the satellites was, untilit was recently dropped, the main cause of errors on the distancemeasurements between the receiver and the satellites. Nevertheless, thecauses of measurement errors remain from the use of the system ofradiopropagation itself, in other words the delays resulting from thepropagation of radio waves through the atmospheric layers, as well asthe errors resulting from possible reflections known as multiple rays.The estimation of these errors and their communication by groundstations, which group two or more measurement receivers and a means ofradiocommunication, allows the user receiving this information tocorrect their own measurements and thus to implement an accuratedifferential positioning. There exist, however, other kinds of errorsthat cannot be eliminated by these differential systems, namely theerrors resulting from a degradation in the operation of the satellitenetwork which can generate measurement errors depending on the physicalcharacteristics of the receivers using these signals. For example, whenthe signal transmitted by a satellite is affected by anomalies known as“EWF3”, the interference alters the correlation mechanisms implementedby the receivers and falsifies the distance measurements in a way thatdepends on the high-frequency filtering analog characteristics of thereceiver and the intervals separating the correlation channels and,consequently, falsifies the determination of the position. Thesephenomena must be detected in order to guarantee the integrity of thedifferential positioning used for landing approach guidance by aircraft.To remedy this, two or more ground stations in the same reception zoneare employed so that, using differential measurements, these errors canbe eliminated or sufficiently attenuated. There exist, however, otherkinds of errors that cannot be eliminated by these differential systems,namely the errors due to the ephemereses providing the satellitepositions and the errors due to the interference affecting the signalstransmitted by the satellites. The errors due to the ephemereses arisefrom the fact that these are manually input by operators and that typingmistakes are always possible. The satellite positioning errors couldthus be several kilometers. In order to eliminate them, it suffices tocompare several consecutive values and eliminate the one that clearlydoes not fit.

When the signal transmitted by a satellite is affected by anomaliesknown as “EWF” (Evil Waveform), the interference distorts thecorrelation peak produced in the ground-based receivers, whichdistortion does not allow the correlation to be performed correctly andtherefore falsifies the determination of position. To remedy this, aworking group has proposed a method for monitoring the quality of asignal transmitted by a satellite in an article numbered WP-13 andentitled “Validation of Revised Signal Quality Monitoring Algorithms forDetecting C/A Code Evil Waveforms” which was presented in Toulouse,France, during the “GBAS Working Group Meeting” which was held from the20^(th) to 24^(th) Mar. 2000 as part of the conference “GlobalNavigation Satellite Systems Panel (GNSSP)”, this working group beingthe “Working Group B”. The method proposed in this article essentiallyconsisted in sampling, at precise points in real time, the correlationfunctions produced in GPS ground-based receivers, in comparing thesesampled values to the set values, and in declaring the received signalinvalid if the result of the comparison exceeded a certain threshold.This method uses precise assumptions about the characteristics of thedetector and the detection criteria are based on the instantaneousobservation of the shape of the correlation peak which entirelydetermines the definition of the receiver and the definition of thedetection algorithms. Another method would consist in systematicallysampling the correlation peak. This method is satisfactory in theory,but, in order to put it into practice, it would require material meansat an exorbitant cost. Indeed, an 18-satellite GPS system, for example,would necessitate 720 correlation channels which would remove anyeconomic advantage of the GPS system, a system which is supposed to beless costly to operate than the existing systems.

SUMMARY OF THE INVENTION

The subject of the present invention is a method for ground-basedmonitoring of the possible presence of anomalies, in particular of theEWF type, in a signal received from a GPS satellite, which method couldbe implemented with the minimum of material means possible at thereceiving station, without however risking the non-detection ofsignificant anomalies in the received signal.

The method of the invention, which is based on the measurement of thedistortion of the correlation peak, consists in taking samples of thecorrelation peak which is produced during the processing of the signalreceived from the satellite, in storing these samples over aninstantaneous sliding time window of at least around 1 minute duration,in storing these instantaneous windows over a period of at least severalhours so as to extract therefrom a statistically determined mean value,in comparing the contents of each instantaneous window to this meanvalue and, if the result is greater than a detection threshold, indeclaring that there is a significant interference affecting thereceived signal and in eliminating the latter.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention will better understood upon reading the detaileddescription of one embodiment, taken as a non-limiting example andillustrated by the appended drawing, in which the single FIGURE is adiagram explaining the weighting step implemented by the method of theinvention.

DETAILED DESCRIPTION OF THE DRAWINGS

The method of the invention applies to a receiver receiving signalstransmitted by geographical positioning satellites, which receiver iscommonly referred to as a GBAS (Ground-Based Augmentation System). Thisterrestrial receiver comprises an SQM (Signal Quality Monitor) functionresponsible for continuously monitoring the quality of the receivedsignals and for warning when the quality is unacceptable in order toreject those signals judged unsuitable for positioning measurements andtherefore to avoid falsifying the measurements.

DETAILED DESCRIPTION OF THE INVENTION

The ground-based station receiver delivers samples of the receivedsignal correlation peak at the rate of two times per second, with a viewto carrying out amplitude measurements in the “in phase” correlator.According to the invention, these samples are five in number and takenat precise instants which are sufficiently characteristic of thecorrelation peak to determine its exact position with the minimumpossible number of samples. These instants are located in a conventionalmanner using relative values with respect to the period of the PNsequence clock frequency, known as “chip”. These values are takensymmetrically with respect to the correlation signal peak, the centralvalue being that of the maximum of the peak (called “prompt”), namely(in values of chip fractions): prompt, ±0.044, ±0.088.

In addition, both short-term and long-term statistical analyses arecarried out on the correlation peaks originating from the signalsreceived from each of the satellites concerned, for each of theaforementioned five values, in order to obtain the individualstatistical characteristics of these values as a function of theconditions of reception of these signals at the ground-based receivingstation. Five standard deviation values σ_(i) (namely: σ_(−0.088),σ_(−0.044), σ_(prompt), σ_(0.044) and σ_(0.088)) and five mean valuesμ_(i) (namely: μ_(−0.088), μ_(−0.044), μ_(prompt), μ_(0.044) andμ_(0.088)) are thus calculated for each type of analysis (short-term andlong-term) and for each satellite concerned.

The short-term analyses are effected within a sliding time window of atleast around one minute duration, and the long-term analyses within asliding time window of at least several hours duration, cumulating allthe short-term analyses relating to this long-term window. A variableΔμ_(i) is then defined, such that Δ_(i)=μ_(i)(ct)−μ_(i)(λt) withμ_(i)(ct) being the value of μ for the sample of rank i considered fromthe peak relative to the short-term analysis and μ_(i)(λt) the value ofthe same sample relative to the long-term analysis. A weighted criterionsqm for the quality of reception at each analysis period is thencalculated, the criterion being given by the following relation:${sqm} = \frac{\begin{matrix}{{\left\lbrack \frac{\Delta\;\mu_{- 0.088}}{s_{- 0.088}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{- 0.044}}{s_{- 0.044}} \right\rbrack^{2} +}} \\{{\left\lbrack \frac{\Delta\;\mu_{prompt}}{s_{prompt}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{+ 0.044}}{s_{+ 0.044}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{+ 0.088}}{s_{+ 0.088}} \right\rbrack^{2}}}\end{matrix}}{MDE}$

In this relation, MDE is a detection threshold analytically determinedso as to obtain a desired false alarm probability ratio (for example,7.2×10⁻⁸ for the OACI standard). If the value of sqm thus calculated isgreater than 1, the presence of an abnormal waveform, or EWF, isdeclared and, consequently, the signals received from the correspondingsatellite must be rejected.

It will be noted that the term (sqm)² follows a chi-square statisticallaw with four degrees of freedom.

With reference to the OACI standard, the equivalent of the Kffdcoefficient, which has a value of 5.26 according to this standard for adistribution with a false alarm probability of 7.2×10⁻⁸, has a value of5.36, in the case of the invention, for a statistical distributionfollowing the chi-square law. Accordingly, owing to the fact that thesqm criterion is weighted, MDE must have a value of 5.36 in order toobtain the same false alarm probability ratio.

A shift register 1 which receives, at one end, the stream 2 of PRN codesof the signal received from a satellite is shown schematically in thesingle FIGURE. The stream of internal PRN codes generated in thereceiver at the ground-based reception station is indicated by an arrow3.

These internal codes have the values that the samples of the correlationpeak should have at the aforementioned sampling instants (central point,±0.044, ±0.088) if the received signals were not affected by parasiticEWF. The internal codes corresponding to instants −0.088, −0.044,prompt, +0.044 and +0.088 are each sent to an input of a convoluter, 4to 8 respectively, whose other input respectively receives the followingvalues: contents of the register for the instant −0.088, contents forthe instant −0.044, difference of the contents of the registers relativeto the instants +0.044 and −0.044 (obtained by a subtractor 10), andcontents for the instants +0.044 and +0.088. In addition, the contentsof the register for the instant where the prompt should appear are sentto a convoluter 9.

The six resulting correlation channels at the output of the convoluters4 to 9 are respectively: I_(−0.088) and Q_(−0.088), I_(−0.044) andQ_(−0.044), I_(Δ) and Q_(Δ) (“delta” mode), I_(0.044) and Q_(0.044),I_(0.088) and Q_(0.088), I_(prompt) and Q_(prompt) (“point” mode). The“delta” and “point” mode channels are used to follow the correspondingsatellite, and the four other channels are used for the I and Q (inphase and in quadrature) measurements at the four corresponding samplingpoints of the correlation peak.

1. A method for the ground-based monitoring of EWF-type anomalies in apositioning satellite signal, characterized in that it consists intaking samples of the correlation peak which is produced during theprocessing of the received signal, in storing these samples over aninstantaneous sliding time window of at least around 1 minute duration,in storing these instantaneous windows over a period of at least severalhours so as to extract therefrom a statistically determined mean value,in comparing the contents of each instantaneous window to this meanvalue and, if the result is greater than a detection threshold, indeclaring that there is a significant interference affecting thereceived signal and in eliminating the latter.
 2. The method as claimedin claim 1, characterized in that the samples are five in number foreach correlation peak analyzed, and in that instants corresponding tothe appearance of the maximum of the peak, to ±0.044 chip and to ±0.088chip relative to the peak maximum are sampled.
 3. The method as claimedin claim 2, characterized in that the weighted criterion (sqm), whosepurpose is to determine the quality of the received signal, is given bythe relation: ${sqm} = \frac{\begin{matrix}{{\left\lbrack \frac{\Delta\;\mu_{- 0.088}}{s_{- 0.088}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{- 0.044}}{s_{- 0.044}} \right\rbrack^{2} +}} \\{{\left\lbrack \frac{\Delta\;\mu_{prompt}}{s_{prompt}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{+ 0.044}}{s_{+ 0.044}} \right\rbrack^{2} + \left\lbrack \frac{\Delta\;\mu_{+ 0.088}}{s_{+ 0.088}} \right\rbrack^{2}}}\end{matrix}}{MDE}$ in which Δμ_(i)=μ_(i)(ct)−μ_(i)(ct), where μ_(i)(ct)is the value of μ for the sample of rank i considered from the peakrelative to the short-term analysis and μ_(i)(λt) is the value of thesame sample relative to the long-term analysis, σ_(i) is the standarddeviation for each of these samples and MDE is a detection thresholdwhich is determined so as to obtain a desired false alarm probabilityratio, and where the value of sqm thus calculated must exceed unity forthe presence of a significant interference to be declared.